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Abstract. Recent anomalies found in cosmological datasets such as the low multipoles of the Cosmic 
Microwave Background or the low redshift amplitude and growth of clustering measured by e.g., abun¬ 
dance of galaxy clusters and redshift space distortions in galaxy surveys, have motivated explorations 
of models beyond standard ACDM. Of particular interest are models where general relativity (GR) 
is modified on large cosmological scales. Here we consider deviations from ACDM-I-CR within the 
context of Horndeski gravity, which is the most general theory of gravity with second derivatives in the 
equations of motion. We adopt a parametrization in which the four additional Horndeski functions 
of time ai{t) are proportional to the cosmological density of dark energy fluEit). Constraints on this 
extended parameter space using a suite of state-of-the art cosmological observations are presented 
for the first time. Although the theory is able to accommodate the low multipoles of the Cosmic 
Microwave Background and the low amplitude of fluctuations from redshift space distortions, we find 
no significant tension with ACDM-I-GR when performing a global fit to recent cosmological data and 
thus there is no evidence against AGDM-I-GR from an analysis of the value of the Bayesian evidence 
ratio of the modified gravity models with respect to ACDM, despite introducing extra parameters. 
The posterior distribution of these extra parameters that we derive return strong constraints on any 
possible deviations from ACDM-I-GR in the context of Horndeski gravity. We illustrate how our 
results can be applied to a more general frameworks of modified gravity models. 
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1 Introduction 

The latest observations of the cosmic microwave background (CMB) by the Planck satellite [1-4] and 
of large structure clustering by several surveys e.g., [5, 6] have provided us with a better understanding 
of the Universe. Indeed, despite the large increase in data volume and accuracy, the ACDM model 
seems to be extremely successful at describing all these data. The ACDM model is based on six free 
parameters that are adjusted to the observations, and assumes homogeneity, isotropy and that the law 
of gravity is Einstein’s general relativity (GR). However, as data have improved, some “inconsisten¬ 
cies”, “tensions”, “anomalies”, (renamed “curiosities” in [1]) with the ACDM model have appeared 
and while they seem to be not statistical significant, they are “persistent”. These “anomalies” are: 
i) the low i multipoles in the observed CMB temperature anisotropy power spectrum are lower than 
predicted by the standard model, ii) CMB (Planck) direct measurements of the lensing potential are 
mildly in tension with the prediction of the angular power spectrum best-fitting ACDM model , Hi) 
the redshift space distortion of galaxy clustering data, the BOSS Lyman—a power spectrum and the 
iv) CFHTLens constraint on the present-day amplitude of fluctuations at linear scales of 8 Mpc/h, 
as, are also lower than predicted by the CMB-calibrated ACDM model, v) The observed cluster 
abundance seems also to be lower than predicted by the model, which could also be explained by a 
lower as- However the level of tension with ACDM depends somewhat on how the cluster sample is 
selected, with Planck SZ-selected clusters giving the stronger tension [7] . Finally vi) the value of 
the Hubble constant inferred from the CMB within a ACDM model is lower than the value directly 
measured locally [8] and e.g., Ref. [9] and refs therein.^ 

These anomalies are each at the « 2cr level and therefore taken individually are not of enough 
statistical significance as to claim any “new physics” beyond the standard ACDM model. However, 
if all these anomalies could be fit by a single model, they would become more than “curiosities” and 
the possible discovery window for new physics opens up. Indeed, there has been already significant 
activity in the literature to try to explain these anomalies 

One basic assumption of the ACDM model is its reliance on Einstein’s general relativity theory. 
One obvious step is then to investigate if changing this assumption could provide a better fit to 
the current observational data. The models proposed in the literature that deviate from Einstein’s 
gravity are abundant, and in principle infinite. Then, in the absence of strong theoretical priors that 
can identify a particular model as the “real” theory of gravity, a generic framework is highly preferred. 
However, increasing the generality of this framework increases also its complexity. It is then crucial 
to balance our theoretical beliefs in order to optimize the freedom that is left to fit the observations. 
To achieve this description many attempts have been done in the literature, such as the parametrized 
post-Friedmann [13] or an Effective Field Theory approach [14, 15]. In the present paper we assume 
that modifications of gravity on cosmological scales are well described by an additional scalar degree 

^The (Tg tension is what drives claims in the literature for massive neutrinos e.g., [7, 10] or phantom dark energy 
[11], the Ho tension drives claims for sterile neutrinos [12]. 

search the ADS database reports in excess of 100 articles addressing this issue. 
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of freedom with at most second-order derivatives in the covariant equations of motion. An additional 
requirement is that the theory should satisfy the weak equivalence principle, i.e. all matter species are 
coupled minimally and universally to the metric This leads to the Horndeski lagrangian [16-18]. 
It encompasses many of the classical dark energy (DE) and Modified Gravity (MG) models studied 
to explain the late-time cosmic acceleration: quintessence, kinetic gravity braiding, galileons, f{R)- ^ 

An efficient description of the linear perturbation theory in the Horndeski lagrangian has been 
investigated in two independent ways in [14, 15, 26-28] and in [29] (see also [30] for an equivalent 
approach at second-order). We will make use of the notation of the latter, but the two results 
are equivalent. The maximal freedom one can have is described by four functions of time besides 
the Hubble parameter H{t), which is responsible for the expansion history of the universe, and one 
constant, i.e., the amount of matter density today. One of the advantages of this approach is that we 
can separate the background evolution from the perturbations. Indeed, these four functions appear 
only at the perturbative level, they are independent from each other and with respect to H(t). Here 
we choose a particular parametrization for these functions proposed by [29] and used in [30], where 
they are proportional to the cosmological density of dark energy This reduces the effective 

freedom of the theory to four free parameters. 

The aim of this paper is to constrain the free parameters of our description using the most recent 
CMB and large scale structure data. 

Horndeski gravity includes GR as a special case, thus if Horndeski gravity provides a signifi¬ 
cantly better statistical fit to the data away from the GR limit this will signal possible new physics. 
Conversely if not such signal is found the results can be used to place limits of possible deviations 
from GR on cosmological scales. 

Within the ACDM model and its popular extensions it is quite widespread to perform global 
fits considering compilations of state-of-the-art cosmological datasets. Specific models, such as / (i?) 
or galileons, have been considered or specihc data sets e.g. [31-40]. This however has not been 
attempted systematically for generic deviations from GR such as Horndeski models. Recently, the 
Planck collaboration [41] studied the implications of some DE/MG models with the latest Planck 
data. In particular, they first parametrized the background evolution only, letting the perturbations 
evolve with standard equations on top of a non-ACDM background. Then, they parametrized directly 
the perturbations introducing two independent functions of the curvature perturbations [42], using 
the public available code MGCamb [43, 44]. Besides the analysis of some “real” theory of gravity, 
such as massive gravity, / (i?) or coupled quintessence, they also studied DE/MG with the EFT 
approach described in the previous paragraph. They used the public code EFTCamb [34, 45], to get 
constraints on a sub-class of the Horndeski class of models, i.e. in our notation the alphas different 
from zero are uk and ae = —om (see [29] and Sec. 2 for details). Their conclusion is that there are 
no significant deviations from ACDM for all the models investigated. This can be ascribed to the 
lack of more precise/accurate data, or to the theoretical priors they are imposing, i.e. not enough 
freedom in the parameter space. While for the first limit we have to wait for the next generation of 
surveys, the latter can be relaxed to cover more generic frameworks. This is the scope of this work. In 
addition here for the hrst time we present a global fit for Horndeski gravity considering a compilation 
of state-of-the-art cosmological data including the latest large-scale structure data. 

The rest of this paper is organized as follows. In §2 we review the Horndeski parametrization 
and outline the methodology adopted and the data sets used. Results and their interpretation are 
reported in §3. Finally we conclude in §4. 

2 Methodology 

We begin by reviewing the Horndeski model and the role of the free functions of the model. Since with 
cosmological data it is not possible to constrain non-parametrically these functions we parametrize 

®Note that with the Horndeski lagrangian it is possible to describe particular classes of models that break the weak 
equivalence principle, such as the ones studied in [19, 20]. However, general frameworks as non-universal disformal 
couplings, or Lorentz-invariance-violating models (e.g., Horava-Lifshitz gravity), or beyond Horndeski theories, e.g., [21- 
25] are not described by the Horndeski lagrangian. These models therefore are not considered here. 


- 2 - 



them and motivate our choice. We also present the data-sets we use here. We consider CMB data 
and large-scale structure data obtained from galaxy surveys. In this paper we do not consider clusters 
data, weak lensing data or local measurements of the Hubble parameter. The amount of tension 
between the ACDM model and cluster abundance depends on the cluster data set chosen and more 
specifically on the mass-observable relation, indicating that this probe is not as mature and robust 
as other probes of We therefore leave the cluster abundance for future work. Traditional -2D” 
weak lensing analysis uses non-linear scales which could be affected by poorly understood physical 
processes (e.g., baryonic effects). Restricting the analysis to linear scales increases significantly the 
error-bars making all “tension” disappear [46-48] . Moreover the minimal description of the Horndeski 
class of models as it has been studied in the literature is strictly valid for linear and second-order 
perturbations [30]. Even if the goal of modeling non-linear scales in generic theories beyond GR is 
in principle achievable, i.e., it is possible to find a finite set of alphas at every order in perturbation 
theory, it is not developed enough and therefore using non-linear scales is well beyond the scope of 
this paper. Finally the tension with the direct Hq measurement turns out not to be too significant 
after recent reanalysis see e.g., [49] and refs therein for discussion. 

2.1 The Horndeski model and adopted parametrization 

The Horndeski action is the most general action for a single scalar field that has second-order equations 
of motion on any background and satisfies the weak equivalence principle. By using an Effective Field 
Theory approach, which has been developed for Horndeski models in [14, 15, 24, 26-29, 50] it is possible 
to identify the minimum number of operators that fully specify the linear evolution of cosmological 
perturbations. The result [29] is that the maximal freedom of the Horndeski lagrangian at linear- 
order in perturbation theory can be described efficiently by five functions of time plus one constant, 
i.e., the amount of matter density today = 0). One of the functions is the Hubble parameter 

H{t), which is responsible for the expansion history of the universe. The other four functions we use 
appear only at the perturbative level, and they were introduced first in [29]. In general Horndeski 
theories ilm{z = 0) can be computed only through measurements of the large scale structure [51, 52]. 
This is due to the fact that at the background level DE/MG can have a component that mimics 
some matter species. As an example, there can be a model that predicts that the DE/MG density 
scales as poEif) = Pdbo + PDEiO,{t)~^, where pdeo and pdei are two free parameters. Then, there 
would be no clear separation between matter and DE just observing the expansion history of the 
universe, and fixing H{t) would not be enough to determine the matter content today. However, this 
degeneracy could be removed at the perturbative level by using large scale structure measurements. 
In this paper, for simplicity we assume that both H{t) and the matter content today are specified by 
a fiducial ACDM model, i.e., the time dependence of ^Ide{z) has no components oc ^^(z). 

Here, it is useful to summarize the physical properties of these functions, in order to understand 
how to parametrize them and to interpret the results of this paper: 

• kineticity: ok. It is the most standard kinetic term present in perfect fluid structures (it is 
the only function different from zero in simple DE models as quintessence). It is possible to 
demonstrate that, if the quasi-static (QS) limit is valid, it does not enter the equations of motion. 
This indicates that it can not be constrained for models with speed of sound close to the speed 
of light [53] . Indeed, in this case the sound horizon of the theory coincides with the cosmological 
horizon and the QS approximation would be validated on the scales of interest. 

• braiding: It comes from a mixing between the kinetic terms of the metric and the scalar. 

The more classical models that have a braiding different from zero are /(i?)/Brans-Dicke and 
cubic galileons. It modifies the growth of perturbations, i.e. / and Geff. As shown in [29], it 
is responsible for all the new scale dependence of the theory beyond the usual Jeans length. If 
ub = 0, it is still possible to have anisotropic stress between the curvature perturbations, but 
Geff — 1. For this reason, this property is crucial in order to modify the shape of the power 
spectrum, and the low i multipoles in the observed GMB spectrum. 

^As studied in [54], the braiding in General Horndeski theories has interesting features at the non-linear level. 
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• Effective Planck-mass and Planck-mass run rate: , aM- A fixed redefinition of the Planck 

mass can not have a detectable effect. Indeed, it can be hidden in the field equations by 
appropriately redefining the densities. For this reason in this work we will not pay too much 
attention to M^. What we can observe is a variation of the Planck-mass, i.e. the Planck-mass 
run-rate. An evidence of om 7 ^ 0, indicates that the theory of gravity is non-minimally coupled. 
As shown in [55], Om generates anisotropic stress in the curvature perturbations and modifies 
the evolution of the gravitational waves. In the presence of non-zero braiding it contributes to 
modify the growth of perturbations through a modification in the effective Newton’s constant. 

• tensor speed excess: a^. It is defined as the excess in the speed of tensors with respect to the 
speed of light. Even if this is its main property, which in principle could be observed [56-58], 
on the scalar sector it introduces anisotropic stress and it modifies the growth of perturbations 
if ae 7 ^ 0 . 


Hence, ob is responsible for all the new scale dependence of the theory beyond the Jeans length. 
If Ob = 0, the effective Newton’s constant, which regulates the relation between the curvature per¬ 
turbation (the trace of the spatial metric) and the matter overdensities, assumes the standard value 
on sub-horizon scales a?H^), see Eqs. (4.9) and (4.11) of [29]. Then, the evolution of linear 

perturbations on sub-horizon scales in the extreme QS limit (which is a usual approximation assumed 
when semi-analytical calculations are performed) can be considered as the evolution of the perturba¬ 
tions in standard perfect fluid structures, with the possibility to have a non-zero anisotropic stress in 
the curvature perturbations. As noticed in [29], it is then possible to define the braiding scale ksit) 
as 


k%{t) 

a{tyH{tY 


-ftmit) + 2 ^ 2 ”*" q ,2 ^ — arit)) . 


( 2 . 1 ) 


Note that the formulation of the braiding scale here is different from the formula given in [29], 
since we are assuming a ACDM background. The ratio between ax and Og regulates the magnitude 
of kg{t). In the absence of the second term of Eq. (2.1), the braiding scale lies at the cosmological 
horizon. But in the presence of non-minimal coupling (i.e. au 7 ^ ot 7 ^ 0), it is sufficient to tune 
Ob ctK in order to push the braiding scales to very small scales (smaller than the scales we are 
observing). At this point one could ask what the magnitude of kg{t) is in a pure ACDM universe. 
Indeed, if we consider a DE model like quintessence (ob = 0), as the closest model to ACDM, then 
—>■ 00 . This would indicate that standard gravity is recovered on super-braiding scales, where 
the equations of motion have a perfect fluid structure. On the contrary, we could reduce the Horndeski 
class of models to a simple MG theory, such as f{R) (ok = 0), and consider it as the closest model to 
ACDM. In the limit where the scalar field associated acquires a large mass it would be decoupled from 
gravity {f{R) —)■ R), leading to standard results. In this example we could conclude that standard 
gravity is recovered on sub-braiding scales. Then, the correct answer to the previous question is 
simply that the braiding scale defines the transition between two different physics. The exact ACDM 
limit has an indeterminate braiding scale. 

When considering standard results in the literature, some confusion arises because, in general, 
non-minimally coupled models (i.e., om 7 ^ 0 or ax 7 ^ 0) have also a non-zero braiding. Then, it 
becomes natural to associate ok with DE, and ob with MG. However, in our framework, all these 
operators are independent from each other, which ultimately means that MG can lie equivalently 
on super-braiding or on sub-braiding scales. We conclude that it is easier to have standard gravity 
on super-braiding scales because it is sufficient to choose au = ot = 0. On sub-braiding scales 
the Planck-mass run-rate and the tensor speed excess have to be tuned to suppress the effect of the 
braiding on the growth of perturbations. 

The description we use has five functions of time. One of them, i.e. H(t), fully specifies the 
expansion history of the universe, while the other four modify the evolution of the linear perturbations. 

However, current data are not sufficient to constrain non-parametrically all these functions of 
time. Therefore we assume ACDM as our fiducial model for the evolution of the background. At 
the perturbative level for deviations from ACDM, we follow the parametrization proposed in [29, 30], 
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Ok 

Ob 

Om & Ot 

(k-Q/aHY 

Gravity type ^ 

^ 0 

^ 0 

^ 0 

Indeterminate 

Quasi-ACDM (GR) 

0(1) 

—>■ oo 

“Strong” MG (non-GR) 

0(1) 

^ 0 

~0(1) 

Minimal MG (non-GR) 

0(1) 

~0(1) 

Classical MG (non-GR) 

0(1) 

^ 0 

^ 0 

—>■ oo 

Classical DE (GR) 

0(1) 

—>■ oo 

Generalized DE (non-GR) 

0(1) 

^ 0 

~0(1) 

Minimal MG (non-GR) 

0(1) 

~0(1) 

Horndeski (non-GR) 


ok/ob 

(om — Ot) /ob 

{ k - Q/aHY 

Gravity type ^ 

^ 0 

^ 0 

~0(1) 

Minimal MG 

0(1) 

~0(1) 

Classical MG 

—>■ oo 

Indeterminate 

“Strong” MG 

—> 00 

^ 0 

Indeterminate 

Quasi-ACDM 

0(1) 

—>■ oo 

Classical DE 

—>■ oo 

—>■ oo 

Generalized DE 



Parameters 

Gravity type 

k'^ > k'/^ 

{ob, om, Q:t} 

MG 

k'^ « k'i 

{ok,om, ax} 

Generalized DE 


Table 1. We show the general behavior of the Horndeski class of models varying the functions Ui or a 
combination of them. In the top table we assume each Oi to be “—>■ 0” or “O (1)”, we compute the estimated 
magnitude of ks and we give the expected type of gravity in the deep sub-horizon limit. In the central table 
we show the magnitude of fee and the expected gravity type as a function of particular combination of the Oi 
functions; this particular combination is closely related to the definition of the braiding scale, Eq. (2.1). In 
the bottom table we give the general behavior and the list of the important parameters on sub/super-braiding 
scales. This clearly indicates that we expect two different physics in these two regimes. 


which reads 

ai(t) = VLde {t) Ci. ( 2 . 2 ) 

Here, all the a functions are fixed to be proportional to the dark energy density thus reducing the 
extra freedom of the theory to four parameters, i.e. {ck, cb, cm, ct}. With this choice we ensure 
the standard evolution of the universe during the eras dominated by radiation and matter, while we 
allow for deviations from standard gravity at late-times. This assumption, even if it is not unique, is 
motivated as follows: The a-functions defined in terms of the Horndeski functions (see Appendix A 
of [29]) are driven by the same functions of the scalar field as the energy density and its derivatives. 
Thus, one would naively expect the parameters Ci to be 0(1). In addition, this parametrization 
describes exactly the behavior of simple shift-symmetric models as the imperfect fluid on its attractor 
[59], and approximately more complicated models as the best fit of covariant galileons [39, 60]. 

In Table 1 we list the various limits of the Horndeski class of models varying the alpha functions. 
In addition, we give a prediction for the magnitude of the braiding scale when it can be estimated. 
The ACDM limit is recovered exactly by setting all Ci = 0. However, ck represents the contribution 
given by perfect fluid structures, i.e., what it is usually classified as DE. Then, even if we set cb = 
Cm = Ct = 0 and leave ck free to vary, the theory of gravity would still be GR. 

2.2 Analysis method 

We use a Bayesian approach to infer the constraints on our model given the datasets described in 
Sec. 2.3. The standard fitting of models defined in a multi-dimensional parameter space is done 
through Monte Carlo sampling techniques. These allow us to infer the posterior of the model’s 
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parameters, which we obtain by sampling the likelihood of the parameters values given the data and 
assuming a uniform prior in those parameters. 

In order to do this, we use the code HiClass® [61] to solve the Boltzmann equations. This 
code, used for the first time in [60] , is an extension of the Cosmic Linear Anisotropy Solver Software 
(CLASS)® by [62] to include Horndeski class of modihed gravity theories. Since as outlined above 
we choose the parametrization in which the time evolution of the four free functions of this class of 
theories is given by the evolution of the constraints quoted in this paper are actually on the 

coefficients that fix the proportionality between a and ^de- 

This Boltzmann code is then interfaced by the cosmology sampler MontePython ' by [63] 
which is a powerful engine to run Monte Carlo Markov chains (MCMC). We run this code using the 
Metropolis-Bastings sampler to obtain our chains until we consider them converged. In this paper 
our convergence criterion is given by a value of i? — 1 ~ 0.03 of the Gelman and Rubin parameter 
[64]. In order to avoid instabilities, we also assume a hard prior on the coefficient ct > —0.9® and a 
initial Planck mass of Mf > 1. We also assume in our runs two massless and one massive neutrino 
with rrii, = 0.06eV. Finally, we remove a small fraction of the points which correspond to a negative 
braiding scale squared = 0). 

2.3 Datasets 

We use a compilation of recent cosmological data including CMB, large scale structure, redshift space 
distortions and baryon acoustic oscillations. To test the robustness of the results, different data sets 
combinations are analyzed, before proceeding to the global fit. 

CMB: We include the Cosmic Microwave Background (CMB) temperature power spectrum from Planck. 
We compare three cases: we consider 1) Planck 2013 temperature power spectrum data [1] with 
the low-£ multipoles of the polarization power spectrum from WMAP [67] (hereafter WP), 2) 
Planck 2015 temperature and polarization data [3, 4], and 3) Planck 2015 temperature and po¬ 
larization data with WMAP polarization replacing Planck polarization data at low multipoles 
(Planck2015WP). At the end we present our main results for our “gold” data set Planck2015WP 
and refer the reader to the Appendix. We choose the latter as our default case because we find 
WP to be more constraining than Planck low-£ polarization for certain non-standard models 
(see Section 3). We emphasize that we use the full shape of the power spectrum rather than 
the shift parameters that determine the position of the acoustic peaks. This dataset includes 
about 7 e-folds out of the ^ 60 that we expect from inflation. This is the only probe we use 
from the high-z Universe, much earlier than when late-time acceleration of the Universe took 
place. On the other hand, the multipoles below £ < 90 probe super-horizon scales, which are 
relevant to this work, especially due to the anomalies reported in these scales (e.g., low value 
for the quadrupole [68, 69], low value around £ ~ 20 [70], etc.). 


P{k): We use the shape of the power spectrum of galaxies from the WiggleZ survey [6] including 
scales in the range 10“‘^h/Mpc < k < 0.2/i/Mpc. These probe roughly another 7 e-folds. This 
constrains the shape of the power spectrum in the late Universe 0.2 < z < 1.0 so this becomes 
sensitive to late-time effects. There are also claims of anomalies with these kind of datasets 
(e.g., [10]), mostly related with the value of the normalization of matter fluctuations erg which 
is slightly different from the one derived from CMB data. Effects from massive neutrinos have 
also been invoked to bring these two datasets into better agreement. Uncertainties from mod¬ 
eling of the halo occupation distribution that maps the relation between the galaxy and matter 
fluctuations have also been claimed [71]. 


®http://hiclass-code.net 
®http://class-code.net 
^http://montepython.net/ 

®For values ct < —0.9 instabilities in the scalar sector prevent our MCMC to explore this region of parameter space. 
This region is anyway disfavoured because observations of binary pulsars place a strong limit on the lower bound of the 
speed of tensors, i.e. cx > —10“^ [65, 66]. 
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Survey 

Measurement 

Value 

Reference 

6dFGS 

rd/Dv{z = 0.106) 

0.327±0.015 

[72] 

SDSS-MGS 

Dv{z = 0.l5)/rd 

4.47±0.16 

[82] 

BOSS-LOWZ 

Dv{z = 0.32)/rd 

8.47±0.17 

[74] 

BOSS-CMASS 

Dv{z = 0.57)/rd 

13.77±0.13 

[73] 

BOSS-LYA 

Da{z = 2M)/rd 

10.8±0.4 

[83] 

BOSS-LYA 

c/H{z = 2.36)rd 

9.0±0.3 

[83] 


Table 2. List of the BAO measurements used in this paper. 


BAO: Another cosmological probe we use is Baryon Acoustic Oscillations (BAO). BAO have their 
origin in the sound waves of the baryon-photon plasma in the early Universe, which were later 
imprinted in the distribution of galaxies, which we have observed with galaxy surveys such as 
6dFGS [72], BOSS [73, 74], or WiggleZ [6]. These show a feature at a characteristic scale of 
^ 150 Mpc that is used to constrain the distance-redshift relation, effectively being used as a 
standard ruler. We note that these measurements tend to marginalize over the shape of the 
power spectrum (or correlation function) and focus on the position of the BAO peak, so they 
cannot constrain by themselves the matter power spectrum. We note that some anomalies have 
been reported regarding the BAO in the distribution of the Lyman-a forest, which seems to 
be in tension with the distance derived from the best-fit ACDM model derived from Planck 
data. Since this BAO measurement corresponds to a redshift of z = 2.4 in which we should not 
expect any influence from late-time acceleration effects, it is interesting to see if they affect our 
constraints on MG. BAO measurements help fix the expansion history thus reducing possible 
parameter degeneracies. The BAO measurements used are summarized in Tab. 2 and shown in 
the left panel of Figure 1. 


RSD: Finally we also include data from Redshift Space Distortions (RSD) in the galaxy distribution. 
These are also derived from Large Scale Structure surveys such as 6dFGS [75], SDSS [76, 77], 
BOSS [78], WiggleZ [79], or VIPERS [80]. Assuming linear theory, it is possible to use the 
anisotropies in the galaxy clustering at intermediate scales {30h~^ — 150/i“^Mpc) to infer the 
line-of-sight contamination in the galaxy redshifts due to peculiar velocities in bulk motion [81]. 
Similarly to (isotropic) galaxy clustering being able to constrain the combination bag, where 
b is the bias of the galaxy distribution with respect to the dark matter and erg the ampli¬ 
tude of (linear) perturbations on 8Mpc//i scales, RSD measures the combination fag, where 
/ = d log D/d log a is known as the growth rate and is the logarithmic derivative of the growth 
factor. It is therefore a probe of the growth of structures and hence of MG. Since these measure¬ 
ments have been obtained with galaxy surveys below redshift z <1, these should be affected by 
the same effects that caused the late-time acceleration of the Universe. Moreover, there have 
been claims that /ug values from RSD are systematically low compared with those obtained 
from the best-fit ACDM model derived from Planck data [4], which makes this dataset very 
relevant for this study. The redshift space distortions data compilation we adopt is summarized 
in Tab. 3 and shown in the right panel of Figure 1. 


3 Results 

We now present our constraints on the Horndeski parameters derived from Monte Carlo sampling of 
this parameter space using the cosmological datasets described in the previous Section. For CMB data 
we use the compilation Planck2015WP. We have also produced constraints using Planck 2013-1-WP 
and using Planck 2015 data [3, 4], which uses Planck’s own polarization measurements instead of 
the low-£ polarization data from WMAP. Because of the longer integration time and the additional 
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Survey 


fiz)<T8{z) 

Reference 

6dFGS 

0.067 

0.423±0.055 

[75] 

MGS 

0.15 

0.53±0.I9 

[77] 

LRG 

0.30 

0.49±0.09 

[76] 

WIGGLEZ 

0.22 

0.42±0.07 

[79] 

WIGGLEZ 

0.41 

0.45±0.04 

[79] 

WIGGLEZ 

0.60 

0.43±0.04 

[79] 

WIGGLEZ 

0.78 

0.38±0.04 

[79] 

BOSS 

0.57 

0.450±0.0II 

[78] 

VIPERS 

0.80 

0.47±0.08 

[80] 


Table 3. List of the RSD measurements used in this paper. 





T3 







redshift redshift 


Figure 1. Left panel: BAO measurements of the angle-averaged distance Dv{z)lrd used in this paper. Right 
panel: RSD measurements of the growth rate times the amplitude of fluctuations f{z)as{z) used in this paper. 
We also show the predictions from ACDM (black solid line) and MG best-fitting models to CMB only (blue 
lines: ck = 0; red lines: ck = 1; green lines: ck = 10). 


information from the high i polarization offered by Planck 2015, the naive expectation is that Planck 
2015 should be more constraining than Planck 2013. If we enforce a hard prior on the integrated 
optical depth due to reionization Treio > 0.04 that is indeed the case. However, our results depend 
on the choice of the prior. We conclude that for the non-standard models considered here the low i 
polarization measurements from WMAP are more constraining that Planck’s. This is why we choose 
to report our updated constraints using as the CMB dataset a combination of the high multipoles of 
Planck temperature and polarization power spectrum, the low multipoles of Planck temperature power 
spectrum, and the low multipoles of WMAP polarization power spectrum. A detailed comparison 
between constraints from these different CMB datasets is presented in the Appendix. 

First of all, we note that the coefficient ck is largely unconstrained by the data, even when the 
sampling is done in logarithmic space. This was expected, since ok does not introduce new scale 
dependence and it does not modify the linear growth. In the quasi-static approximation where time 
derivatives are considered to be sub-dominant with respect to space derivatives, ok does not enter 
the equations of motion. In this regime we expect no effect of ok on the observables and thus no 
sensitivity of our analysis to this parameter. Therefore, we choose a few arbitrary fixed values of ck 
and report constraints for these cases. In Table 4 we report the 95% confidence level bounds on the 
values of cb,m,t using different combinations of datasets and after marginalizing over all the other 
parameters in the MCMC chain. In interpreting this table one should bear in mind that the case 
Ck = cb = Cm = ct = 0 corresponds to the ACDM model. 

The table shows that there is some mild dependence of the constraints on the adopted value of 
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_cb_Cm_ 

CMB, ck = 0 +0.04 < Cb < +1.91 —0.86 < cm < +2.00 

CMB, CK = 1 -0.01 < CB < +1.84 -0.77 < cm < +2.00 

CMB, CK = 10 +0.14 < CB < +1.92 -0.73 < cm < +2.00 

CMB+BAO, CK = 0 +0.07 < cb < +1.95 -0.83 < cm < +2.00 

CMB+BAO, CK = 1 +0.08 < cb < +1.96 -0.85 < cm < +2.00 

CMB+BAO, CK = 10 +0.19 < cb < +1.97 -0.79 < cm < +2.00 

CMB+RSD, CK = 0 +0.27 < cb < +2.48 -1.43 < cm < -0.23 

CMB+RSD, CK = 1 +0.15 < cb < +2.37 -1.39 < cm < -0.13 

CMB+RSD, CK = 10 _ +0.20 < cb < +2.39 -1.40 < cm < -0-14 

CMB+PK, Ck = 0 +0.00 < cb < +1.82 —0.85 < cm < +2.00 

CMB+PK, Ck = 1 —0.03 < cb < +1.85 —0.83 < cm < +2.00 

CMB+PK, CK = 10 +0.12 < CB < +1.87 -0.80 < cm < +2.00 

CMB+BAO+RSD+PK, ck = 0 +0.24 < cb < +2.32 -1.36 < cm < -0.13 

CMB+BAO+RSD+PK, ck = 1 +0.10 < cb < +2.29 -1.35 < cm < -0.08 

CMB+BAO+RSD+PK, ck = 10 +0.19 < cb < +2.30 -1.36 < cm < -0.06 


_Ct_ 

-0.90 < CT < +1.20 
-0.90 < Ct < +1.08 
-0.90 < Ct < +0.77 
-0.90 < CT < +1.14 
-0.90 < CT < +1.21 
-0.90 < CT < +0.87 
-0.90 < CT < -0.44 
-0.90 < CT < -0.44 
-0.90 < CT < -0.44 
-0.90 < CT < +0.99 
-0.90 < CT < +0.88 
-0.90 < CT < +0.72 
-0.90 < Ct < -0.39 
-0.90 < Ct < -0.41 
-0.90 < Ct < -0.41 


Table 4. Constraints on the coefficients cb, cm, and ct from different cosmological dataset combinations and 
for different values of ck. Quoted limits are 95% CL. A hard prior on ct > —0.9 is applied as well as a prior 
on —2 < Cm < +2 that has become relevant in some cases. 


Ck. In all cases (all data set combinations and all values of ck) the parameter cm is consistent with 
zero. For this parameter the ck dependence is very mild. 

When RSD is not included, the ct parameter is consistent with zero and shows no ck depen¬ 
dence for the data set combinations CMB, CMB+BAO and CMB+PK. For the CMB+RSD and 
CMB+BAO+PK+RSD data sets this parameter is negative and 0 is formally outside the 95% con¬ 
fidence. The Ct lower limit is imposed by the prior ct > —0.9. It is then clear that RSD pushes ct 
towards negative values. Indeed, of the data we consider RSD is the only one that probes (with good 
statistical significance) the growth of perturbations (in fact it measures /erg). There is some small 
dependence of that in CMB leasing and ISW but at lower significance. It is well known that within 
the LCDM paradigm, RSD wants an amplitude of fluctuations slightly lower than that predicted when 
fitting the model to CMB data. Negative values of the parameter ct achieve just that: a lower /erg 
at the redshift of interest for RSD without major alterations to other observables. 

The confidence region for the parameter cb also excludes 0 in few cases and positive values seems 
to be formally favoured. The upper limit of cb does not depend on the adopted value of ck, but it is 
sensitive to the choice of data set. The lower limit on the other hand is also sensitive to the adopted 
Ck value. 

Our tightest constraints, which include all the datasets considered, formally appear to favour a 
theory of gravity that is minimally coupled with tensors propagating slower than the speed of light. 
The results are robust with respect to changes in the low-£ CMB polarization data (from Planck to 
WMAP, see Appendix 4). We have also tested that the model-dependent RSD measurement in [78] is 
not driving our results, by checking the consistency of the constraints when that data point is replaced 
with the measurement in [84]. These tests confirm the robustness of our results. 

For the other data sets the ACDM model (ck = cb = cm = ct = 0) is not always included in 
the 95% posterior confidence interval. 

We will study below the Bayesian evidence ratio of the ACDM to MG models. But before we do 
so, in order to interpret the significance of this result we examine carefully the values of the effective 
statistic -given by the log of the likelihood- measured at the best-fit point in our Monte Carlo 
chains. This is shown in Table 5, where we also include the standard ACDM case (in which the MG 
equations are not used) for direct comparison. We can see that CMB data actually prefers a larger 
value of the coefficient ck, and with a lower significance this is also true for RSD data. On the other 
hand, PK data slightly compensate this trend by selecting lower values of ck at the best-fit model. 
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- log likelihood 


CMB 

(low-.«) 

CMB 

(high-^) 

BAO 

RSD 

PK 

Total 

A Total 
^ACDM 

CMB (ACDM) 


1014.7 

1222.2 




2236.9 

- 

CMB (MG), CK = 0 


1014.4 

1221.9 




2236.3 

0.6 

CMB (MG), CK = 1 


1014.3 

1221.8 




2236.1 

0.8 

CMB (MG), CK = 10 


1013.5 

1221.4 




2234.9 

2.0 

CMB-tBAO (ACDM) 


1015.1 

1221.6 

5.6 



2242.3 

- 

CMB-tBAO (MG), ck = 0 


1014.8 

1221.8 

5.0 



2241.5 

0.8 

CMB-tBAO (MG), ck = 1 


1013.9 

1221.1 

5.2 



2240.2 

2.1 

CMB-tBAO (MG), ck = 10 


1014.0 

1221.7 

4.8 



2240.5 

1.8 

CMB-tRSD (ACDM) 


1014.3 

1223.7 


4.6 


2242.6 

- 

CMB-tRSD (MG), ck = 0 


1013.7 

1222.1 


2.2 


2237.9 

4.7 

CMB-tRSD (MG), ck = 1 


1013.9 

1223.3 


2.1 


2239.3 

3.3 

CMB-tRSD (MG), ck = 10 


1013.4 

1222.3 


2.7 


2238.4 

4.2 

CMB-tPK (ACDM) 


1015.3 

1221.5 



228.6 

2465.5 

- 

CMB-tPK (MG), CK = 0 


1014.1 

1222.6 



229.8 

2466.5 

-1.0 

CMB-tPK (MG), CK = 1 


1013.8 

1222.4 



229.2 

2465.4 

0.1 

CMB-tPK (MG), CK = 10 


1013.3 

1221.9 



229.5 

2464.7 

0.8 

CMB-tBAO+RSD-fPK (ACDM) 


1014.5 

1226.5 

4.7 

3.1 

228.2 

2477.0 

- 

CMB-tBAO+RSD-fPK (MG), ck 

= 0 

1013.5 

1222.9 

4.8 

1.9 

229.0 

2472.1 

4.9 

CMB-tBAO+RSD+PK (MG), ck 

= 1 

1013.7 

1222.9 

4.9 

1.8 

228.9 

2472.2 

4.8 

CMB-tBAO+RSD-kPK (MG), ck 

= 10 

1013.7 

1222.7 

4.8 

1.7 

228.9 

2471.8 

5.2 


Table 5. Absolute value of the log likelihoods (i.e. X^/2) at the best ht point from the individual data 
that comprises each dataset combination explored in our analysis. The column labelled Total displays the 
maximum likelihood value in the chain. The last column shows the difference in Log likelihood with respect to 
the ACDM model. Red (negative) numbers represent worst fit, positive (black) numbers better fit. Given the 
intrinsic uncertainty of the MCMC in determining the best likelihood value, the improvement in offered 
by the more complex model is in most cases not significant. 


The MCMC procedure is not optimized to find the best fit model which maximizes the likelihood, 
therefore there is an intrinsic error associated to these numbers which have been estimated to be 
- 0.7[85]. 

Compared to the ACDM model, we find that the improvement on the fitting of cosmological data 
due to the extra degrees of freedom provided by the Horndeski parameters is not significant in most 
of the cases. Possible exception are the CMB (all ck values) and the CMB+BAO for ck = 10 cases 
where the improvement is A log Likelihood ~ 2 at the “cost” of three extra degrees of freedom. This 
suggests that the deviations found in our datasets are still consistent with a fluctuation within the 
ACDM scenario, even though (remarkably) the posterior distributions of MG coefficients presented 
in Table 4 are not always consistent with zero. 

The best-fit models that correspond to the Log Likelihood values of Tab. 5 are shown in Figure 2 
where we present the cases in which we fit only CMB data, CMB and one of the Large Scale Structure 
datasets (BAO, RSD, or PK), and the full combination CMB-I-BAO-I-RSD+PK. We show how the 
observables fitted in our analysis are reproduced by the MG models described in Section 2 where we 
fix the coefficient ck to 0 (blue lines), 1 (red lines), and 10 (green lines). We also include for direct 
comparison the best-fit model in a ACDM cosmology with standard gravity (black lines). We find 
that MG has more flexibility to fit the observables, but, given current error bars, the difference cannot 
be regarded as more than a mere statistical fluctuation. It remains to be checked with future data 
whether the significance of these differences increases or is still consistent with noise. 

Given the above apparent contradiction, we investigate whether the posterior confidence inter¬ 
vals might be driven by strong correlations between the parameters varied in our MGMG chains. 
Strong correlations are known to introduce prior volume effects in multidimensional analysis when 
the posterior distribution is non-Gaussian. 

Interestingly we find no significant correlations between the ACDM and the MG sectors of our 
parameter space, which greatly simplifies this analysis. We will therefore focus on possible cor- 
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Figure 2. CMB temperature power spectra, matter power spectra, and structure growth for ACDM and 
MG best fit models. From top to bottom, we show these quantities when we ht only the CMB, CMB+BAO, 
CMB+RSD, CMB+PK, and CMB+BAO+RSD+PK respectively. Vertical dashed lines correspond to the 
braiding scale fee (a = 0) for each particular model and dataset combination. 


relations between the parameters of the MG sector. In figure 3 we show the 68% (light shade) 
and 95% confidence level contours (dark shade) between pairs of MG parameters when the full 
CMB+BAO+RSD+PK dataset combination is used. This figure shows that, especially in the cb 
case, it is the long non-Gaussian tail at high values that drives the effect. The constraints on the MG 
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Figure 3. Constraints on cb, cm, ct and M* from the combination of CMB+BAO+RSD+PK datasets. 


parameters are nevertheless found to be uncorrelated, at least to the extent driven by the precision 
of the datasets used here. 

The issue of whether the data favours a MG model is an issue of model selection or model 
comparison rather than parameter estimation. Within the Bayesian framework, which is the one 
adopted here, model comparison is done by considering the model-averaged likelihood, referred to 
here as the Bayesian Evidence, E. Under the assumption of equal a priori model probabilities, the 
ratio of Evidence values for two models, given the same data, quantifies the relative odds of these 
models being the correct description of the observations. Thus the key model comparison quantity 
is the Bayes factor, which is the ratio of the Evidence values for the two models. In general, the 
Evidence is the result of a multi-dimensional integral over the model parameters, the evaluation of 
which can be computationally expensive. However here the two models are nested: the simpler model 
(the ACDM-I-GR one) is a specihc case of the MG models when the four coefficients cb,m,t,k are zero. 
In this case it is possible to perform rigorous model comparison between nested models without the 
need for evaluating numerically-intensive multi-dimensional integrals using the Savage-Dickey density 
ratio (SDDR; [86, 87]). 

In Table 6 we report the Bayes factor of the AGDM to MG models computed following [85]; we 
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Dataset combination 

Evidence ratio In {Ej^cdm/Eb) 

interpretation 

CMB 

1.48 

substantial 

CMB-l-BAO 

-0.21 

not significant 

CMB-hRSD 

0.64 

not significant 

CMB-hPK 

1.13 

substantial 

CMB-hBAO-hRSD-hPK 

0.09 

not significant 


Table 6. The Savage-Dickey Density Ratio for the ACDM -|- GR model with respect to the Modified Gravity 
models studied here. We have considered the case ck = 0. 


use a slightly modified version of the Jeffrey’s scale to interpret the evidence ratios. The Bayes factor 
favours the simpler, ACDM, model in the cases CMB and CMB-hPK with odds ^3:1 (“substantial” 
evidence). 

As with ah Bayesian methods, the Bayes factor between the two models depends on the prior on 
the model parameters. The comparison between the Bayesian Evidence ratio and the Log Likelihood 
analysis presented above, which is prior-independent, serves to quantify a possible prior-dependence 
of the Bayes factor. 

Next we briefly discuss the implications of our bounds on the Ci parameters on the existence 
of a (possible) new scale associated to gravity (eq. 2.1). While formally one can compute a value 
for ksiz = 0 ) given a set of parameters Ci, the interpretation of confidence levels on this parameter 
requires some caution. As shown in Table 1, the ACDM limit (which is recovered exactly by setting ah 
our parameters Ci = 0) has an indeterminate fee. In our analysis the ACDM model offers a very good 
fit to the data. For the collection of models in parameter space close to ACDM, the distribution 
is unconstrained. When the ACDM model is not disfavoured by the data the braiding scale is not 
a good (derived) parameter to use since small variations of our Ci can produce huge changes on this 
scale. For this reason we do not dwell much on this quantity here. The braiding scale can still be a 
good quantity to use when investigating models beyond the simple ACDM model: it is meaningful to 
evaluate it only in a MG scenario. 

At this point it is useful to give a prescription on how to interpret our results for “real” theories. 
Indeed, strictly speaking the constraints we obtain are only valid for sub-classes of Horndeski that 
have the time evolution of their a’s proportional to the time evolution of the DE/MG density on 
a ACDM background. In the literature, it is possible to find many models that are sub-classes of 
the Horndeski lagrangian, but in general each of them has a different dynamics. We claim that our 
results can be applied to more general frameworks. The reason for this is that most of the DE/MG 
models are constructed in such a way that they explain the late-time acceleration of the universe, with 
standard evolution at early-times (radiation/matter dominated eras) and modifications at late-times. 
Then, one can have in a reasonable amount of time (without the need of running MGMG chains) an 
idea of what is the allowed/not allowed region in the parameter space of a given theory. As a toy 
model we will make use of a minimal generalization of the imperfect fluid introduced in [59, 88 ]. In 
this model we consider K = — X”, G3 = G4 = Mpj/2 and G 5 = 0.® This is a shift-symmetric 

model that has an attractor for the background evolution. As usual, on this attractor the equations 
can be integrated analytically, which is an advantage we use in order to illustrate some results. The 
steps one has to follow are: 

1. Solve the background evolution (i.e. H{t)) and check that it is close to AGDM. Within our 
example it is possible to demonstrate that in the infinite future, where ^de = 1 , the equation 
of state for the scalar field is wde = — 1 ; 

2. Use the definition of the alphas in terms of the Horndeski functions (in [29]) in order to express 
them as a function of background quantities. In our case we get: ax = 6 n(l — 2n -|- 2m)flDE, 
a-B = 2nVt£iE and ctm = ax = 0 ; 

®These functions are the usual Horndeski functions appeared first in [17]. The exact notation we used in order to 
define them can be found in Sec. 2, precisely Eq. (2.1), of [29]. The model parameters therefore are n and m. 
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cb Cm ct 

CMB+BAO+RSD+PK, ck = 0 -0.01 < cb < +3.06 -1.65 < cm < +0.22 -0.90 < ct < -0.13 

CMB+BAO+RSD+PK, ck = 1 -0.12 < cb < +2.91 -1.61 < cm < +0.27 -0.90 < ct < -0.14 

CMB+BAO+RSD+PK, ck = 10 +0.00 < cb < +2.86 -1.59 < cm < +0.30 -0.90 < ct < -0.15 


Table 7. Constraints on the coefficients cb, cm, and ct from different cosmological dataset combinations and 
for different values of ck. Quoted limits are 99.73% CL. A hard prior on ct > —0.9 is applied. 




Cb Cg c^ 


Figure 4. Correlations between the coefficients Cb, Cm, and ct- Contours are shown for 68.3%, 95.4%, and 
99.7% CL. The symbols correspond to different values for the parameters of the model of Ref [59, 88] as 
discussed in the text. The three points in the left two panels correspond to (from left to right) n = 1/3, 
n = 2/3, n = 1. In the right panel the three models overlap. Clearly values n > 2/3 are disfavoured. 


3. Check that the evolution of these a functions is indeed (X^Ide (or in which regime this applies). 
In our case this condition is satisfied exactly; 

4. Invert the relations between our Ci and the parameters of the theory (usually this step has to 
be done numerically). In our the result is: n = cb/ 2 and m = (ck/cb + cb — l)/2; 

5. Compare with constraints in Fig. 4 to define the allowed parameter region of the original theory. 

In Fig. 4 we illustrate this last step in the procedure. The symbols correspond to different values 
for the parameters of the toy model considered here. Values of n ^ 1/3 or 2/3 are allowed but values 
n > 3/4 are excluded. This is interesting as n = 1, which is the standard kinetic term and standard 
case proposed in the original reference, is ruled out. 

4 Discussion and conclusions 

We have considered a general version of Horndeski theories of gravity, where modifications to general 
relativity on cosmological scales are well described by an additional scalar degree of freedom with at 
most second-order derivatives in the equations of motion and where the weak equivalence principle 
is satisfied. This description encompasses many of the classical DE/MG models studied to explain 
the late-time cosmic acceleration. We chose a parametrization in which the time evolution of the four 
arbitrary functions in the Horndeski Lagrangian is given by the time evolution of ^de- The ACDM 
limit is recovered when all the extra parameters are zero. 

We have performed a global fit for the parameters of the model (including the cosmological param¬ 
eters), considering the state-of-the art cosmological data: Cosmic Microwave Background data from 

“It is important to note that this condition can be satisfied approximately even in more general frameworks as the 
best fit models of the covariant galileons [39, 60]. 
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the Planck mission and the matter power spectrum, redshift-space distortions and Baryon Acoustic 
Oscillation data from galaxy surveys of Large Scale Structure. 

We used the standard Bayesian approach to parameter inference and obtained posterior confi¬ 
dence intervals on the model’s parameters. The main result of this work is that we find no significant 
statistical evidence for deviations from Einstein’s gravity in the current data. While formally for some 
data-sets combinations the ACDM-I-GR limit is (just) outside the 95% posterior confidence interval, 
the improvement in the fit at the expense of adding extra parameters, quantifies in terms of difference 
of log likelihood is not significant. The Bayes factor also confirms this by not favouring the more 
complex model (Horndeski gravity) over the simpler (ACDM-I-GR) one. These findings are robust to 
to choice of specific data-set used and deviations from the ACDM-I-GR limit are mostly driven by the 
low multipoles of the Cosmic Microwave Background anisotropies. 

These results are the best constraints so far on this parametrization of the Horndeski Lagrangian, 
in which ai oc which limit any possible deviations from ACDM-I-GR. 

It is interesting to note that in constraining this class of models beyond ACDM, the measurement 
of large-scale Cosmic Microwave Background polarization plays an important role. In fact we find 
that the low £ polarization data from the Planck satellite are less constraining than the WMAP 
ones. Forthcoming improvements on large-scale polarization measurements might help improve our 
constraints. 

While strictly speaking the constraints we obtain are only valid for sub-classes of Horndeski 
considered, our results can be applied to a more general framework. In fact most of the DE/MG models 
reported in the literature are constructed in such a way to explain the late-time acceleration of the 
universe, with standard evolution at early-times (radiation/matter dominated eras) and modifications 
at late-times, making therefore possible to “map” them to our class of model by interpreting our 
parameters as effective parameters. This mapping is of course model-dependent and might not always 
be analytic, but this procedure enables one to rule in or out specific models or parameter ranges for 
models without having to perform any data analysis or expensive exploration of multidimensional 
parameter-space and inference. We have illustrated this procedure for a specific MG model. 
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-0.13 

-0.90 

< 

Ct 

< 

-0.39 

Planckl3-kWP-kBAO-l-RSD-kPK, ck = 

= 1 


-t 0.02 

< 

Cb 

< 

-t2.75 

-1.52 

< 

Cm 

< 

-kO.29 

-0.90 

< 

Ct 

< 

- 0.21 

Planckl5-klowTEB-kBAO-kRSD-kPK, 

CK = 

1 

-tO.34 

< 

Cb 

< 

-t2.90 

-1.57 

< 

Cm 

< 

-k0.04 

-0.90 

< 

Ct 

< 

-0.25 

Planckl5-|-tauprior-|-BAO-fRSD-|-PK, 

Ck = 

1 

-0.03 

< 

Cb 

< 

-t 2.01 

-0.91 

< 

Cm 

< 

-1-0.37 

-0.90 

< 

Ct 

< 

- 0.01 

Planckl5-klowT-tWP-|-BAO-tRSD-|-PK, ck 

= 1 

-to.io 

< 

Cb 

< 

-t2.29 

-1.35 

< 

Cm 

< 

-0.08 

-0.90 

< 

Ct 

< 

-0.41 

Planckl3-kWP-kBAO-kRSD-kPK, ck = 

= 10 


-t0.13 

< 

Cb 

< 

-t2.75 

-1.51 

< 

Cm 

< 

-1-0.35 

-0.90 

< 

Ct 

< 

- 0.22 

Planckl5-klowTEB-l-BAO-kRSD-kPK, 

CK = 

10 

-tO.35 

< 

Cb 

< 

-t2.94 

-1.58 

< 

Cm 

< 

-kO.Ol 

-0.90 

< 

Ct 

< 

-0.24 

Planckl5-|-tauprior-|-BAO-|-RSD-|-PK, 

Ck = 

10 

- 0.02 

< 

Cb 

< 

-tl.97 

-0.89 

< 

Cm 

< 

-kO.36 

-0.90 

< 

Ct 

< 

- 0.02 

Planckl5-klowT-tWP-kBAO-tRSD+PK, ck 

= 10 

-t0.19 

< 

Cb 

< 

-t2.30 

-1.36 

< 

Cm 

< 

-0.06 

-0.90 

< 

Ct 

< 

-0.41 


Table 8. Constraints on the coefficients Cb, Cm, and ct from different Cosmic Microwave Background datasets 
in combination with BAO+RSD+PK and for different values of ck- Quoted limits are 95% CL. A hard prior 
on Ct > —0.9 is imposed. 


Appendix A: Constraining power of different CMB datasets 

In this Appendix we justify our default choice for the CMB dataset combination. In Table 8 we show 
how our constraints change when we choose a different CMB dataset combination. Although the 
Planck 2015 case seems inconsistent with the others (especially for cb), it also forces the optical depth 
to reionization r^eio to be less than 0.04, which we know to be inconsistent with other cosmological 
observations such as e.g., the Gunn-Peterson effect (see e.g., [89]) or WMAP data themselves. We con¬ 
clude that the most constraining CMB dataset while consistent with other cosmological observations 
is the one used in our default case. 
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